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Abstract 



It has recently been suggested that cosmologically significant numbers of 
black holes could form during a first-order QCD phase transition. Further, 
it has been asserted that these black holes would have masses corresponding 
naturally to the inferred mass (~ l-M©) of the MACHOs responsible for the 
observed gravitational microlensing events. In this model, the underlying 
spectrum of primordial density perturbations provides the fluctuations that 
give rise to black holes at the epoch of the QCD transition. We employ 
a simplified model to estimate the reduction in the critical overdensity of 
a horizon-sized primordial perturbation required for collapse to a black hole. 
We find that a first-order QCD transition does indeed produce a sharp peak in 
the black hole mass spectrum, but that this peak corresponds to the horizon 
mass at an epoch somewhat earlier than the cosmological transition itself. 
Assuming a COBE normalized primordial density perturbation spectrum with 
constant spectral index, for the black holes so produced to be cosmologically 
significant would require an extremely finely tuned "blue" primordial density 
perturbation spectrum. Specifically, in the context of our simplified model, 
a spectral index in the range n = 1.37 — 1.42 corresponds to the range 0, ~ 
10 -5 — 10 3 of the black hole contribution to the present-day density parameter. 

97.60.Lf, 12.38.Aw, 95.35.+d, 98.35.Gi 



1 



Typeset using REVT^X 



I. INTRODUCTION 



Recent studies of gravitational microlensing suggest that about half of the Galaxy's halo 
mass can be accounted for by massive compact halo objects (MACHOs) with mass ~ O.5M [] 
HJ. However, a direct search for halo dwarfs in the Hubble Deep Field has limited the halo 
mass fraction of red dwarfs brighter than magnitude Mi = 15 to /rd < 0.06, and a halo 
mass fraction of white dwarfs with Mj < 15 to /wd < 0.33 ||. Comparing what a present 
day white dwarf Galactic halo fraction should have looked like in the past with deep galaxy 
surveys leads to a limit of /wd < 0.10 [|J. A study of the infrared background and Galactic 
metallicity shows that it is likely that /wd < 0.25 [Q. 

While white dwarf halo mass fractions which exceed these limits might be obtained with 
specially tailored assumptions [§-0] j ultimately it may be necessary to consider more exotic 
candidates for the MACHOs. In particular, Jedamzik has recently proposed that horizon- 
size black hole formation from primordial density fluctuations could operate particularly 
efficiently during a first-order QCD phase transition ||. A first-order confinement transition 
involves the coexistence of quark-gluon and hadron phases |§ . On large scales, the pressure 
response of this "mixed" phase is greatly reduced: when compressed, the energy density does 
not increase from increases in pressure and temperature, which remain constant. Instead, 
low density hadron phase is exchanged for high energy density quark-gluon phase. Because 
of the reduced pressure response, the required overdensity for black hole formation is smaller. 
For fluctuations entering the horizon during this epoch, statistically more abundant, lower 
amplitude perturbations could form black holes, leading to a peak in the primordial black 
hole (PBH) mass spectrum at the QCD epoch horizon mass. Since the horizon mass is 

M h (T) « O.87M ( ) 2 ( -^-) ~ V2 , (1) 

v ; ° VIOOMeV/ V51.25/ ' V ; 

where g* is the effective number of relativistic degrees of freedom at cosmological temperature 
T, PBHs formed at the QCD transition temperature Tqcd ~ 100 MeV could be in the right 
mass range to be MACHOs. This possibility takes on added interest in view of the claim 
that galactic halos filled with ~ O.5M black holes could bring about enough nearby BH- 
BH mergers to cause a few events/year in the first generation network of gravitational wave 



interferometers 10 



Other than a general guess that the critical overdensity for collapse would be reduced 
by a factor of order unity during the QCD transition PJllfj, no quantitative estimate of this 



reduction factor has appeared. A relativistic numerical calculation is in preparation [|12 
The purpose of this paper is to use a simplified model to estimate this reduction factor, and 
to assess the cosmological significance of the resulting peak in the PBH mass spectrum. 



1 The absence of short timescale microlensing events restricts the MACHOs to masses ^ 0.1M@, 
unless thick disk/light halo models of the Milky Way are employed. This eliminates brown dwarfs 
unless, in conjunction with these nonstandard Galactic models, the brown dwarf mass function 
peaks close to the hydrogen burning limit p]. It has also been suggested that the gravitational 
lenses are not in the Galactic halo at all, but rather in "the warped and flaring" Galactic disk [||. 
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II. CONDITIONS FOR PBH FORMATION 



A. Radiation equation of state 



A heuristic criterion for gravitational collapse to a black hole of an overdense region of 
energy density p and size S was given by Carr and Hawking These authors argued 

that collapse would ensue when the "gravitational energy" ~ Gp 2 S 5 was large enough to 
overcome the "kinetic energy" of expansion and the pressure forces (here G is the gravita- 
tional constant). The overdense region expands more slowly than the average cosmological 
expansion, and eventually turns around and begins to collapse. Carr and Hawking consider 
the region at the moment of turnaround (when the kinetic energy of expansion is zero), 
and account for the pressure forces by comparing the gravitational energy to the "internal 
energy." This "internal energy" is approximated as 



where p is the pressure and an equation of state of the form p = wp is assumed (w is 
a constant). Then the condition for an overdense region to overcome pressure forces and 
collapse to a black hole is Gp c S 2 > w, where the subscript c indicates that the quantity is 
to be evaluated at the moment of turnaround when the collapse "begins." This condition 
can also be written S c ^ (w/Gpc) 1 ^ 2 ~ Rj, where Rj is the Jeans length. However, the 
overdense region cannot be larger than the horizon size at turnaround, or it will constitute 
a separate universe, causally disconnected from our own, and would not be relevant to this 
discussion. 

For radiation (w = 1/3), the size of the overdense region at turnaround, its Schwarzschild 
radius, the Jeans length, and the cosmological horizon size are all of the same order of 
magnitude. Therefore it is expected that an overdense region of radiation satisfying Gp c S 2 ^ 
1/3 will form a black hole of order the horizon mass, with little chance for centrifugal or 
turbulent forces to have a significant effect. For dust (w=0), the Jeans length is much smaller 
than the horizon size, and centrifugal and turbulent forces could prevent full collapse to a 
black hole. If a hole does form in this case it could be much smaller than the horizon mass. 

The condition for an overdense region to collapse to a black hole is expressed above in 
terms of the density and size of the region at turnaround. However, cosmological density 
perturbations are typically classified by specifying their overdensity 5 = (p — p)/p when 
they enter the horizon, where p is the energy density of the overdense region and p is the 
average cosmological density. For assessing the cosmological significance of the black hole 
population, it will therefore be convenient to express the condition for collapse in terms of 
5, which in this paper will always refer to the overdensity at horizon crossing. 

To make this connection we follow Carr ||14|| . The line element for the "average" universe 
can be taken to be the conformally flat Friedmann-Robertson- Walker (FRW) form, 



pS 3 ~ wpS 3 , 



(2) 



ds 2 = -dt 2 + R 2 [dr 2 + r 2 (d6 2 + sin 2 Od(f)% 



(3) 



where the evolution of the time- dependent cosmic scale factor R is given by 




(4) 
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Here the "average" background geometry/evolution is obtained by viewing the universe on 
a large enough scale that the (nearly critical) mass-energy density can be approximated as 
homogeneously and isotropically distributed. 

We consider "top-hat" fluctuations: spherical, homogeneous overdense regions. Such 
overdense regions, if taken to have homogeneously distributed internal mass-energy, can be 
described by the metric of a closed FRW universe, 

ds 2 = -dr 2 + S 2 [(l - nr^dr 2 + r 2 (d6 2 + sin 2 6d<J)% (5) 

with the evolution of S given by 

8nG 




pS 2 - k, (6) 

where p is the energy density in the overdense region. Obviously, at some point in the 
gravitational collapse of such a fluctuation, pressure gradients will build up and it will not 
be possible to approximate the internal mass-energy distribution as homogeneous. However, 
this approximation will not be bad until near turnaround. Here S can be interpreted as a sort 
of "radius" of the evolving fluctuation. Likewise, r is a timelike coordinate characterizing 
the evolution of the overdense region. 

The general relativistic gauge freedom inherent in this model is fixed by setting = t^, 
Sh = Rh, and (dS/dr)h = (dR/dt) h , ensuring that the when the overdense region enters 
the horizon it is only a density perturbation. (The subscript h denotes a quantity evaluated 
when the overdense region under consideration enters the horizon.) At horizon crossing, the 
energy density of the fluctuation is ph = ph(l + 5). Using these conditions, and employing 
an equation of state p = wp, the evolution of the overdense region is given by 




1 + 5 5 



S 1+3w Rtf 



3ui 



(7) 



so at turnaround 



sl+3w = Rl+3w ( i+£\ _ (8) 



5 

The quantity appearing in the collapse condition is 

G Pc S 2 c = Gph(1 ^f h(1+W) = ( ) 
where we have used Eq. ([8p. Using Eq. (|J) and the fact that the horizon size ~ t, 

^-©r© 2 - 1 ' (10) 

so that the condition for collapse is 

5>w, (11) 

or 5 ^ 1/3 for radiation. In other words, in this conventional picture, the minimum fluctua- 
tion amplitude 5p/p at horizon crossing required for collapse to a black hole is of order unity, 
i. e. quite large. This large required initial amplitude simply reflects the "hard" equation of 
state of radiation. 
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B. First-order QCD transition 



It is not at all clear at this point whether or not there is a first order cosmic vacuum phase 
transition associated with the QCD epoch. The chiral symmetry breaking event remains the 
prime candidate for such a transition, though the lattice gauge calculation results on this 
issue remain unsettled. For illustrative purposes, we will assume that there is a first order 
phase transition associated with this epoch in the early universe. The QCD era is promising 
for primordial black hole production, since phase separation at this time implies that the 
duration of a first order phase transition will be (essentially) a gravitational timescale. 

A bag-like picture employing a QCD vacuum energy provides a simple description of a 
first order QCD phase transition with critical transition temperature Tqcd- The pressure 
and energy density of the hadron phase (T < Tqcd) are approximately 

1 7T 2 

Phad(T) = -phad, Phad(T) = —g had T 4 , (12) 
with, for example, #had ~ 17.25. The quark-gluon phase (T > Tqcd) is described by 

1 7T 2 7T 2 

Pqg( T ) = 3^^s T4 _ 5 ' = ^9 qg T 4 + B, (13) 

where g qg ~ 51.25 and B is the QCD vacuum energy density. During the transition (T « 
Tq CD ), the requirement p ha d(^QCD) = ^(Tqcd) yields 

1 7T 2 

B = 3 30^ qg ~ ^ had ) T Q CD - ( 14 ) 
Averaged over the coexisting phases, the energy density during the transition is 

(p) = /qgPqg(7QCD) + (1 - /qg)Phad(T Q CD), (15) 

where / qg is the volume fraction in quark-gluon phase. The evolution of the average energy 
density during the transition as a function of cosmic scale factor R can be described as 
"dust-like" (p~iT 3 ) §, 

/ R \ 3 / 1 \ 1 

(P){R) = (Pqg + gPhadj - gPhad, (16) 

where R\ is the scale factor when the cosmic energy density reaches p\ = p qg (TQco)- This 
is the start of the phase transition. 

For the purpose of assessing the probability of black hole formation during the QCD 
epoch, we will consider the following further caricature of the behavior of the plasma through 
the transition. For energy densities greater than pi, take the plasma to be radiation with 
9* = #qg- This ignores the vacuum energy density B, which at T = Tqcd contributes about 
18% of the total. For energy densities in the range p2 = Phad(TQco) < P < Pi, we can treat 
the plasma as dust. Here p\ and pi correspond to the energy densities in the plasma at the 
start and conclusion of the phase transition, respectively. Even though the pressure itself is 
not zero, the almost negligible pressure response of the plasma during the transition leads 
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to the dust-like evolution noted above. For densities less than pi y we treat the plasma as 
radiation with = g had . 

The duration of the transition can be obtained by comparing the cosmic scale factor 
at the beginning [R\) and end (R2) of the transtion. From the conservation of co-moving 
entropy density we obtain (R 2 /Ri) = (^qg/ghad) 1 ^ 3 ~ 1-44, so the transition is fairly short. 
An overdense region that goes through the transition during its collapse does not correspond 
to the cases discussed in the previous subsection, i.e. either radiation or dust throughout 
the collapse. How should the condition for collapse to a black hole be modified in this 
case? Based on the above caricature of the behavior of the plasma around the transition, 
we suggest modifying the "internal energy" of Eq. (fj) to (l/3)p c S^ (1 — /), where / is the 
fraction of the overdense region's "evolution volume" spent in the dust-like transition epoch. 

For example, letting Si and S2 respectively denote the size of the initially expanding 
overdense region at the beginning and end of the phase transition, / = (Sf — Sf)/S^ for 
an overdense region that has completed the transition before turnaround. Put another way, 
for an overdense region that obeys p = wp for much of its evolution and a dust-like phase 
for the rest, we define an effective Jeans length Rj )C s = (w/Gp c ) 1//2 (l — f) 1 ^ 2 at turnaround. 
In terms of the overdensity of the fluctuation at horizon crossing, the condition for collapse 
can be expressed 

5>w(l-f). (17) 

Therefore, the larger the fraction of a fluctuation's evolution time which is spent in the mixed 
phase regime, the smaller the required initial horizon crossing perturbation amplitude for 
collapse to a black hole. The evaluation of this collapse criterion will be the subject of Sec. [TV]. 
Note that in our analysis we have assumed that at no time is the universe ever dominated by 
the QCD vacuum energy B. Were such a vacuum-dominated condition to arise, then a kind 
of "mini-inflation" occurs, where the scale factor can increase exponentially with time |T5 



However, in this case the nucleation rate for bubbles of low temperature phase will increase 
very rapidly with time and decreasing temperature, with the result that the fraction / will 
likely be smaller than in the non-vacuum-dominated case. 



III. CLASSIFYING PERTURBATIONS IN THE QCD TRANSITION EPOCH 

"Top-hat" perturbations are completely specified by 5 — their overdensity when they enter 
the horizon — and by the cosmological time at which horizon crossing occurs. Here we will 
find it convenient to use the average cosmological energy density at horizon crossing ph 
rather than the time of horizon crossing. In particular, we use the variable x to identify 
the epoch of horizon crossing, where x = ph/ Pi- Again, p\ is the energy density above 
which the plasma is pure quark-gluon phase, and p2 is the density below which the plasma 
is pure hadron phase. We define the constant f3 = pi/ pi- Then perturbations with x > 1, 
< x < 1, and x < (3~ l enter the horizon when the average cosmological background — 
not the overdense region itself — is in quark-gluon, mixed quark-gluon/hadron, and hadron 
phase, respectively. 

In addition, our description of the collapse of an overdense region involves two important 
events: (1) the region's horizon crossing, and (2) its turnaround, when it stops expanding 
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and begins to recollapse. The state of the matter at these two events constitute various 
classes of perturbations; these are listed in Table |. 

Given the characteristics of an overdense region, 6 and x, we can specify the class A-F 
to which it belongs. The classes of perturbations that exist for various epochs of horizon 
crossing and values of S are shown in Table [H]. From the information in this table two curves 
can be drawn in the (x,5) plane (Fig. [I]). Classes A, B, and C are separated from classes 
D and E by the curve 5 = x~ l — 1, and classes D and E are separated from class F by the 
curve 5 = f3~ l x~ x — 1. These divisions are based on the state of matter in the overdense 
region when it enters the horizon. 

In contrast, the curves in the (x, 5) plane separating classes A and B, B and C, and D and 
E are based on the state of matter in the overdense region at turnaround. For example, the 
boundary between classes A and B is determined by p c = p±, that is, by the condition that 
an overdense region entering the horizon in the quark-gluon phase just reaches the density 
for the phase transition to begin when it begins to recollapse. The condition 

_ p h (l + 6)Rj _ 

Pc - ^4 - Pi, (1°) 

using Eq. (^) with w — 1/3 yields the curve defined by 

<rk='~ l - (19) 

The other curves can be obtained in a similar manner, with the complication that changes 
in the density evolution between radiation (p oc S~ 4 ) and dust (p oc S~ 3 ), and/or vice- versa, 
must be followed when the density reaches pi and/or p 2 respectively. The analysis leading 
to Eq. (H) can then be generalized in a straightforward manner. The curve defining the 
boundary between classes B and C is 

,V> (3- l x~V\ (20) 



(1 + 5) 3 / 2 

and the boundary between classes D and E is defined by 

<5 3 



sy 



(3- l x-\ (21) 



These last two curves are valid only in the regions above and below the curve separating 
classes A, B, and C from classes D and E, respectively. 

Fig. [I] shows the (x, 5) plane divided into regions that constitute classes A-F. Class A 
has no dust-like phase, and 5 > 1/3 for collapse. Class D does not exist for 5 ^ 1.8, so 
perturbations of this class would be expected to collapse even without the dust-like phase. 
Therefore, interesting reductions in the critical value of 5 for collapse to a black hole are 
possible for classes B, C, E, and F. 

For reference, the expressions for the turnaround radius S c for the various classes are 



given in Table |TJ. We also give expressions for S± and S2, the values of S for which the 



density in the overdense region reaches pi and p 2 respectively. These are 
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S 1 = R h x^(l + 5) 1/4 , 
which is valid (and useful) for classes B and C, 

S 2 = R h {{3x) l '\l + 5) l l\ 
valid (and useful) for class E, and 

S 2 = R h (f3x) l '\l + 5) l l\ 

for class F. 



(22) 



(23) 



(24) 



IV. COMPUTING CRITICAL OVERDENSITIES 



In the framework of our model, overdense regions with 



*£g(l-/) 



(25) 



will collapse to a black hole, where / is the fraction of "evolution volume" that the overdense 
region spends in the dust-like phase. For class B, turnaround occurs while the overdense 
region is still in the mixed phase: 



/ = ^ f ^ = l-(|) 3 (Class B), 
and 5*i and S c are given in Eq. (|2|) and Table |T| respectively, yielding 



(26) 



Si 
S c 



3 X 3/2 S 3 



(l + <5) 3 / 2 

For classes C, E, and F, the transition to hadron phase is completed by turnaround, so 

C3 c3 
°2 °i 



(27) 



/ 



55 



(Classes C, E, F). 



For class C, it is convenient to write this in the form 



Mr 

where (S 2 / Sx) 3 = pi/p 2 = [3, and 



Si 



s 2 ^ 3 



i ] 1 



(Class C), 



(28) 



(29) 



x 



3/402 



(Class C). 



.Sj /3V2(1 + 5 )3/4 

For classes E and F it is convenient to write Eq. fl28f ) in the form 



(30) 



s 2 ^ 3 



s 2 



(Classes E, F), 



(31) 
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where (S 1 /S 2 ) 3 = P' 1 and 

'S 2 \ 3 ((3x) 1/2 S 3/2 



S C J (1 + 6) 

S 2 \ 3 ({3x) 3 ^5 3 / 2 



SJ (1 + 5)3/4 



(Class E) (32) 
(Class F), (33) 



using Eqs. (|23|j24]) and Table |TT|. 

For a given x — that is, for overdense regions entering the horizon at a given time — the 
values of 5 for which collapse to a black hole will occur can be found by plotting both 
(1 — /)/3 and S itself as functions of 5, taking care to use the expression for / appropriate 
to the class specified by (x, 5) (see Fig. [1]). Figure |2] shows such a plot for x — 2, from which 
it is clear that collapse will occur for all values of 5 above a certain critical value. However, 
as seen in Fig. || for larger values of x (x = 15 in Fig. [|) an interesting thing occurs: the 
range of 5 for which 6 > (1 — /)/3 has an upper as well as lower bound. Figure f| indicates 
the region in the (x, 5) plane for which collapse to a black hole occurs. The range of x for 
which there is a reduction in the critical overdensity for collapse has an upper bound; and 
while it has no lower bound, the reduction becomes slight at small values of x. 



V. COSMOLOGICAL SIGNIFICANCE 

The potential cosmological significance of primordial black holes lies in the fact that 
the energy density of a black hole population redshifts as -R -3 , in contrast with radiation, 
whose energy density redshifts as R~ 4 . Thus the relative contribution of black holes to the 
cosmological density parameter Q increases as the universe expands. See, for example, the 
work by Crawford and Schramm |16| . 



For black holes formed at a particular epoch, the contribution to the present-day density 
parameter is 



'BH 



2 x 10 7 e(T) ( ( ) h~\ (34) 

V ; V2.75K7 VIOOMeV; ' V ; 



where e(T) is the fraction of radiation energy density converted into black holes at cosmo- 
logical temperature T, To is the present-day microwave background temperature, and h is 
the present epoch Hubble parameter in units of 100 km s _1 Mpc -1 . Assuming that the mass 
of the black holes is the horizon mass at T, we can write 

e (T) = I*" F(S, T) dS, (35) 

where F(S, T) is the probability of a horizon volume at T to have overdensity 5, and S a and 
6b delimit the range of overdensities that contribute to the black hole population. 

The power spectrum (spatial Fourier transform) of density perturbations is typically 
assumed to have a power law form, oc k n , where k is the wavenumber and n is the 
spectral index. Density perturbations generated during inflation typically have n ~ 1, with 
a normal distribution of overdensities: 
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F = 7sW xp ("sfM))' (36) 

where a(M) is the variance at horizon mass M and can be taken to be [ 17 

(l-n)/4 

\ 1O 22 M " 



a(M) = 9.5 x 10- 5 1 ,87) 



when normalized by COBE measurements of the microwave background anisotropy. 

The critical value of 5 above which collapse to a black hole can occur [5 a in Eq. (j35l) 1 is 
given by the curve in Fig. 0J, taking the lowest branch when it becomes multivalued at the 
higher values of x. Noting that perturbations with 5 > 1/3 make negligible contributions to 
Q, we take 5b = 1/3, except for the higher values of x for which the the curve in Fig. (H) is 
multivalued. For this region, 5b is given by the middle branch of the curve in Fig. (f|). 

It is convenient to change from our variable x = Ph/pi, hi terms of which the curve in 
Fig. (f|) was calculated, to the horizon mass M. This connection is made as follows. We 
account for the fact that we have ignored the vacuum energy density B in the quark-gluon 
phase by defining an effective transition temperature TQCD,eff by 



712 -4 D _ *\ 



..gfi'qg^QCD + B — ggPqg^QCD.eff- (^8) 



Using this equation and Eqs. (|l|) and (14), we then have the following relation between x 
and M: 



x 



^qgT QC D,cff V 51 - 25 / V°- 87M o/ VlOOMeV 



) ■ (39) 



For a given value of the spectral index n, e(M) can be computed, where the horizon mass 
M has replaced T in Eq. (^) as the variable indicating the cosmological epoch. 

For the values of n we will have to consider (n ~ 1.4), e(M) is a very strongly peaked 
function, justifying the assumption in Eq. (0) of black holes being formed at a single epoch. 
This sharp peak comes about as follows. For n > 1, there is increasing power on small scales 
[Eq. ([37])]. Also, 5 a decreases with decreasing M (increasing x, see Fig. f|). Both of these 
tend to make e(M) increase rapidly with decreasing M. However, when M gets too small, 
the range of overdensities contributing to black holes is cut off completely [the cusp at x ~ 50 
in Fig. (01)], and e(M) drops abruptly. 



To obtain e for use in Eq. (34), we choose a value of n and find the peak value of e(M). 
The mass of the black holes is the value of M for which e has its peak value. Figures (|5|) and 
(§) show ^bh as a function of n for two choices of transition temperature Tqcd- The fine 



tuning in n required for ^bh ~ 1 is immediately apparent ||11|| . For a given Tqcd the masses 
of the black holes do not vary significantly over the ranges of n displayed: ttibh — 0.11M Q 
for T QCD = 100 MeV, and m BH ^ O.O3M for T QCD = 200 MeV. 



VI. CONCLUSIONS 

Given the potential difficulties in accounting for the observed MACHO population with 
red or white dwarfs, it is interesting to follow up on the suggestion || that the MACHOs 
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could be black holes formed during a first-order QCD confinement transition. We have per- 
formed a semianalytic analysis using a simplified model. This model involves consideration 
of spherical "top-hat" overdense regions and the use of the Friedmann equation to calculate 
their evolution. A caricature of the behavior of the equation of state near the QCD transi- 



tion is employed, along with an extension of the condition given by Carr and Hawking [13 
for primordial overdense regions to collapse to black holes. 

In the context of a COBE normalized primordial density perturbation spectrum charac- 
terized by a power law with constant spectral index n, we find that for the black holes so 
produced to be cosmologically significant would require n ~ 1.4, with extreme fine tuning 
(see Figs. ^-||) PJITfl. While such a spectral index is not in conflict with observations of 



microwave background anisotropies [|18|,|l£||, it conflicts with limits on n based on, for ex- 
ample, evaporation of lower mass black holes formed at earlier epochs (e.g., Refs. |T7|J20|] 
and references therein). In addition, we find that the holes most efficiently formed by this 
mechanism entered the horizon at an epoch somewhat earlier than the time the average 
universe goes through the QCD transition. Accordingly, we find black holes masses in the 
range 0.11 — 0.03M Q for QCD transition temperatures of 100 — 200 MeV. This mass range 
is somewhat low for the estimated O.5M MACHO mass. 

Our model involves some gross simplifications. The condition for collapse to a black 
hole is heuristic. The use of spherical "top-hat" perturbations is not realistic to begin with, 
and using the Friedmann equation for their evolution after they enter the horizon is even 
less realistic. The bag-like model is not an overly detailed description of the behavior of 
matter near the confinement transition, and we have even employed a caricature of the bag 
model. We have assumed, in keeping with the standard lore, that the mass of black holes 
formed during a radiation dominated era will be roughly the horizon mass at the epoch the 
overdense region enters the horizon. 

The results from our simplified model are not particularly encouraging for this mechanism 
of black hole formation, and refinements in the treatment would probably make matters 
worse. Early numerical treatments of black hole collapse had different results, with one set of 
calculations [21] finding black holes to have about the horizon mass, while othere calculations 



22| found that the holes formed with masses somewhat smaller than the horizon mass. 
More recent calculations |23j suggest that, contrary to the usual assumption, the masses 
of primordial black holes formed in a radiation dominated era may be significantly smaller 
than the horizon mass. These same calculations also found that the overdensities necessary 
for collapse in a radiation era are 2-3 times higher than the 8 « 1/3 from the standard 
lore. (Larger values of 5 are also obtained if one keeps various factors of order unity in 
the heuristic analysis leading to the condition 5 > 1/3 for collapse.) Finally, in assessing 
the cosmological significance of the black holes we have assumed a Gaussian distribution 
of overdensities. Non-gaussian effects on the large overdensity tail are likely to be skew- 
negative [TIJ], although they do not weaken the effect on the needed value of the spectral 



index too much [[L7j] 



The shortcomings of our model may make it unsuitble for precise quantitative predic- 
tions, but we believe two conclusions can be drawn from it that would be confirmed with 
refined treatment: 1) The perturbations gaining the most "benefit" from a first-order QCD 
transition enter the horizon somewhat earlier than the time the averaged universe goes 
through the transition. (This effect is also being seen in numerical calculations [|24|].) This 
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means that the mass of the black holes formed by this mechanism will be somewhat smaller 
than the horizon mass when the averaged universe goes through the transition. 2) If density 
perturbations can be characterized as a COBE-normalized power law with constant spectral 
index, a very finely tuned, blue spectrum would be required for the black holes produced by 
this mechanism to be cosmologically significant f8yTT|j. 
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FIGURES 




0.1 1 10 100 

x 

FIG. 1. Regions in the (x, 5) plane corresponding to the classes of perturbations listed in Table 
|. The variable x identifies the epoch a perturbation enters the horizon: x > 1 (x < 1) corresponds 
to overdense regions that enter the horizon before (after) the average cosmological density begins 
the transition from quark-gluon to hadron phase. The quantity 5 is the overdensity 5p/p of the 
perturbation at horizon crossing. 




FIG. 2. Solid curve: (1 — f)/3, for x = 2. Dashed curve: 5. Collapse to a black hole occurs for 
values of S for which the dashed line is above the solid curve. The intersection corresponds to the 
value of 5 at x = 2 in the curve of Fig. El 
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FIG. 4. The curve in the (x, 5) plane indicating which parameter values lead to collapse to a 
black hole. Without the phase transition, this would be a straight line at 5 = 1/3. 
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FIG. 5. The contribution to the present day density parameter by black holes formed during 
the QCD transition, as function of the spectral index of the primordial density fluctuations. The 
transition temperature is taken to be Tqcd = 100 MeV, and the resulting black holes have mass 
» O.llMpj. 




1.36 1.37 1.38 1.39 1.4 1.41 
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FIG. 6. Same as Fig. ||, but for transition temperature Tqcd = 200 MeV. The black holes in 
this case have mass « O.O3M . 
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TABLES 



TABLE I. Classification of over dense regions according to the state of matter at horizon cross- 
ing and turnaround. 



Class Horizon crossing a Turnaround a 

A qg qg 

B qg m 

C qg h 

Dm m 

Em h 

F h h 



a gq: quark-gluon phase, m: mixed phase, i.e. coexistent quark-gluon and hadron phases, h: 
hadron phase. 



TABLE II. Classes of perturbations that exist for various epochs of horizon crossing and ranges 



of 5. 



Epoch of orizon crossing 


Overdensity 


Class 


x > 1 


6 > 


A, B, C 


/r 1 < x < 1 


S > x- 1 - 1 


A, B, C 




< S < x- 1 - 1 


D, E 


x < /r 1 


5 > x- 1 - 1 


A, B 




fi^x- 1 < 5 < x~ l - 1 


D, E 




< 5 < f3^x- 1 


F 



TABLE III. Radius of the overdense region at turnaround. 



Class Radius at turnaround 

a Sc =R h (^y /2 

B S c = R hX -^£±<fl 

C S c = R h ^(lf) 1 ' 2 

D S c = R h £p 

E S c = i4(/3*)V6(I+^ 

F S c = R h (^) 1/2 
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